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Abstract 



In this work we show how to define the action of a scalar field in a such 
a way that Robin boundary condition is implemented dynamically, i.e., as a 
consequence of the stationary action principle. We discuss the quantization 
of that system via functional integration. Using this formalism, we derive 
an expression for the Casimir energy of a massless scalar field under Robin 
boundary conditions on a pair of parallel plates, characterized by constants 
ci and C2. Some special cases are discussed; in particular, we show that for 
some values of ci and C2 the Casimir energy as a function of the distance 
between the plates presents a minimum. We also discuss the renormalization 
at one-loop order of the two-point Green function in the X(f)^ theory submitted 
to Robin boundary condition on a plate. 
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I. INTRODUCTION 



The Casimir force between two uncharged macroscopic bodies in vacuum is widely re- 
garded as arising from the zero-point fluctuations intrinsic to any quantum system. In the 
case of two flat parallel plates a distance a apart the pressure is proportional to a""^ [1]. 
In recent years several groups have performed high-precision measurements of the Casimir 
force between a flat plate and a spherical surface (lens) or a sphere [2], and also between 
two parallel flat plates [3]. In the last case the original Casimir formula was confirmed to 
15% accuracy. 

Due to its fundamental character the Casimir effect has applications in many areas of 
physics, ranging from the theory of elementary particles and interactions [4,5] to atomic and 
molecular physics [6]. Besides, it has analogues in condensed matter physics, for instance 
in fluctuation induced forces [7] and boundary critical phenomena [8]. More recently, its 
relevance to the design and operation of micro- and nano-scale electromechanical devices 
has been emphasized [9]. 

Usually the details of the interaction between the vacuum fluctuations of the quantum 
field and the macroscopic bodies are neglected, and replaced by classical boundary conditions 
(BC) at the boundary of the latter. While Dirichlet and Neumann BC have been extensively 
studied over the past, the more general case of Robin BC has attracted httle attention. 

A field (p is said to obey Robin boundary condition at a surface E if its normal derivative 
at a point on E is proportional to its value there: 

d 

— (f){x) = c(l){x), a; e E. (1) 

Neumann and Dirichlet boundary conditions are particular cases of Robin boundary condi- 
tion: the first one corresponds to c = 0; the other is obtained in the limit c — > oo (assuming 
that dn4> is bounded). 

The mixed case of Dirichlct-Robin (DR) BC were considered in [10] for a 2D massless 
scalar field as a phenomenological model for a penetrable surface, with playing the role 
of the finite penetration depth. Recently, the Casimir energy for a scalar field subject to 
Robin BC on one or two parallel planes was computed in [11]. 

Here we also compute the Casimir energy for a scalar field under Robin BC on two parallel 
planes. However, we introduce a rather different approach which seems more amenable to 
an eventual computation of radiative corrections in models containing interactions. Its 
starting point is the introduction of suitable boundary terms in the action, which allows us 
to compute the partition function of the system without the explicit imposition of Robin 
BC on the fields. In spite of that we show that the two-point Green function does satisfy 
those boundary conditions. We also find agreement with the main results of Ref. [11] for 
the Casimir energy, which was computed there using the more conventional approach of 
summing the zero-point energy of the normal modes of the field. In particular, we show 
that for the mixed case of DR boundary conditions the Casimir energy as a function of a 
develops a minimum, i.e., there is a configuration of stable equilibrium. 

Finally, we study the renormalization at one-loop order of the two-point Green function 
in the A0"^ theory submitted to Robin BC on a plate. Our analysis differs from previous ones 
[8,12] in two aspects: (i) we keep c arbitrary, instead of considering only the the particular 
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cases c = or c — oo, and (ii) we perform the regularization entirely in momentum space. 
This procedure avoids deahng with distributions and test functions, which are unavoidable 
in the mixed coordinate- momentum space regularization used in [8,12]. 



II. THE MODIFIED ACTION 



Let us consider, for simplicity, a real scalar field living in the half-space z > 0,^ satisfying 
the (Euclidean) equation of motion 



and subject to Robin boundary condition at 2; = 0, 



9^0 - c0 



z=0 



= 0. 



(2) 



(3) 



(We shall assume that c > 0, in order to avoid the possible appearance of tachyons in the 
theory.) One can easily verify that Eqs. (2) and (3) are consequences of the stationary action 
principle applied to the Euclidean action 



S[cf>] 



dz 



1 



where 



f0 fL/2 

d'^x := lim lim / dx^ / dx^ 

L-*oo Jo J-L/2 



1 



L/2 
-L/2 



(x,0) 



dx 



Indeed, computing 5S :— S[(j) + rj] — S[(f)] up to second order in 77 we obtain 

SS = I dl'x |(-9,0 + c0) 77 + dz [-d^ + U'{ct>)] Tyj + 0{rf), 



(4) 



(5) 



(6) 



which implies (2) and (3) if is a stationary point of S. 

Until now, we have been talking about a c/assica/ field. What happens when one quantizes 
the theory? In the usual functional integrals approach one has to integrate over all field 
configurations obeying certain boundary conditions. If such is the case, is it necessary 
to retain the surface term in the action? Bordag et al. [13] argue that it is, in order to 
ensure the Hermiticity (more precisely, the self-adjointness) of the fluctuation operator := 
—d^ -\- U"{4>c), where 4>c is the solution to Eqs. (2) and (3). Saharian [14] has also argued 
in favor of such a surface term: without it, the vacuum energy evaluated as the sum of the 
zero-point energy of each normal mode of the field does not agree with the result obtained 
by integrating the vacuum energy density. 

Here we propose a different approach. We retain the surface term in the action, but 
we shall not impose any boundary condition at 2; = on the field configurations to be 



Conventions: h = c = 1, x = {k, z), where x := {x^, . . . ,x'^ ^) and z := x^. 
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integrated over. Somewhat surprisingly, if we treat U{(f)) as a perturbation, the two-point 
Green function of the unperturbed theory does satisfy Robin BC at z — 0, i.e.. 



as^ 



(a, - c) {cPix) 0(xO)o = 0. (7) 
Let us prove it. First of all, we write the partition function of the unperturbed theory 

/ m exp(-^o), (8) 

where 5*0 is given by Eq. (4) without U {</)). Note that we are integrating over all field 
configurations satisfying the boundary condition 0(x, 0) = 0i(x), and then we integrate over 
all configurations of the surface field 0i(x). In other words, we integrate over all possible 
boundary conditions at z = 0. 0(x, z) also satisfies periodic BC in the x-coordinates, with 
period (3 in the a;°-direction and L in the others (note, however, that we have the limits 
/3 — s> oo and L ^ oo in mind). 

We now decompose as a sum of two fields: (f) = (f)Q + rj, where (f)Q satisfies 

a2</.o(x) = 0, (/)o(x,0) =(/)i(x), (/)o(x,oo) =0. (9) 

Note that, because of the boundary conditions imposed on 00, rj vanishes at z — and 
z — oo. Eq. (9) can be solved using Fourier transform; the result is 

^e^■^•>l(k)e-'=^ (10) 

where /c = |k| and <^i(k) is the Fourier transform of 0i(x). 

In terms of the fields (fi and rj, the partition function Zq becomes the product of two 
independent functional integrals: Zq — ZaZb, where 

Za^ /[i^<^i]exp|-/^i(c + A;)(^i(k)(^i(-k)|, (11) 



Zr^ 



[ [Dn]e^p\- [d'x TdzUd^rif]. (12) 
J?)(x,o)=o { J Jo 2 ) 



Let us now compute the two-point Green function. As a consequence of the factorization 
of Zq, one has {(j){x) (f){x'))o = {(i>o{x) (f)o{x'))A + v{^')) b- According to Eq. (11) we have 

{Mk)^rik'))A = ^S^'\k + k'). (13) 

Gombining this result with Eq. (10) we obtain 



^We follow here a procedure very similar to the one employed in Ref. [15] in the context of quantum 
field theory at finite temperature. 
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d'^^ ik-(x-x') ^ 



-k{z+z') 



(27r)^ 



c + k 



On the other hand, 



where satisfies 



(-af + k^) Dr,{k; z, z') = 5{z - z'), D^{k; 0, z') = D^{k; oo, z') = 0. 
One can easily verify that the solution to (16) is given by 

Dn(k;z,z') = — sinh(A;2;<) exp{—kzy), 
where z^{z>) — min(max){2;, z'}. Collecting terms, we finally obtain 



^-k{z+z') 

+ — sinh(A;2;<) exp(— A;2;>) 



c + k k 

One can easily verify that (^(x) 0(a;'))o indeed satisfies Robin BC at z = 0, Eq. (7). 



(14) 



(15) 



(16) 



(17) 



(18) 



III. CASIMIR ENERGY 

Let us now apply our procedure to the computation of the Casimir energy of a free 
massless scalar field subject to Robin BC on two parallel plates located at the planes 
z — Q and z — a, 



z=0 



0, dz(t) + C2^ 







(C1,C2 > 0). 



The Euclidean action for such a system is given by 



I d'^xi dz 



2 



+ ^ci(/)'(x, 0) + ^C2(/)'(x, a] 



and its partition function is given by 



//•</)(x,a)=</>2(x) 
[D(l)i][D(P2] / W] exp{-S). 

J</.(x,0)=</.i(x) 

The Casimir energy can be extracted from Z using the identity 



(19) 



(20) 



(21) 



Eq — — lim — In Z. 



(22) 



As we did in the previous Section, we shall write (p as the sum of two terms: (p = (po + r], 
where 0o is the solution to the classical equation of motion, d'^cpo = 0, that obeys the 



5 



boundary conditions 0o(x, 0) = 0i(x) and 0o(x, a) — 02 (x) (consequently ry(x, 0) = rj{'x, a) — 
0). By Fourier transforming in the x-coordinates one can explicitly solve for (f)Q, obt^^i^i""' 

/)o(x, z) = j -j—^ ^i^^hla f*^^^^-* sinh A;(a - z) + (^2(k) smh. kz] 



J {ZT^}"- smn/ea 

where /c = |k| and <^j(k) is the Fourier transform of 0j(x), j = 1, 2. 

Expressing S[4)] in terms of 0o and 77, we obtain S — Sa-\- Sb-, where 



(23) 



Sa^ j d'^x|/Jd^[^(a^0or] +^Ci02(x,O) + ^C20g(x,a)| (24) 



d X I dz 



(25) 



Since 0o is a functional solely of 0i and 02, the partition function Z can be written as the 
product of two terms: Z — ZaZb, where 



Za^ ip0i]p02] exp(-5^), 



|-T/(x,a)=0 

Zb = / [Dr}] exp{-SB)- 

^r/(x,0)=0 



(26) 



It follows from Eq. (22) that the Casimir energy is given by the sum of two terms, Eq = 
Ea + Eb, of which the second one is the Casimir energy of a field subject to Dirichlet 
boundary conditions on the plates. Since it is a well known result [16], we shall just quote 
the result: 



EB{a) 



■^rf^l (4.)-(^+^)/^c(^+i). 



(27) 



Let us now compute Ea- Inserting (23) into (24) we can rewrite Sa in terms of ipi and 
(P2 as 



(28) 



Mjj(k) = (cj + /c coth ka) 6ij — k csch ka{l — 6ij). 
Changing the variables of integration in Za to (pj, we obtain 

Za = f[Dip,][Dip2] exp{-SA[^i,ip2]) = ndet-^/=^M(k), 



(29) 



(30) 



hence 



EY^{a)^\un -^lndetM(k) 



d'^k 
(27r)« 



In 



C1C2 + A;^ + (ci + C2)/ccoth/ca 



(31) 
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The expression above diverges, hence requires renormahzation. This is achieved by sub- 
tracting from it the quantity 



SEa = hm E\^'^{a) 



d-l 



I 



(2^ 



In[(ci + A;)(c2 + A;)], 



(32) 



which can be interpreted as part of the self-energy of the plates. (An analogous subtraction 
is necessary in the calculation of Eb-) Since SEa does not depend on the distance between 
the plates, it does not contribute to the force between them. Its subtraction from E^^'^^ is 
thus permissible as long as one is interested — as we are — only in the Casimir force. The 
result of the subtraction is given by 



EA{a) 



-d-l 



I 



{2ny 



In 



1 + 



2(ci + C2)k 1 
{ci + k){c2 + k) ' e^^"- - 1 



(33) 



Performing the angular integration and adding the result to Eb, Eq. (27), we finally obtain 
the Casimir energy for the massless scalar field under Robin boundary conditions: 



-Eo(ci,C2;a) = EB{a) + 



-d-l 



(47r)<^/2r(i) 



dkk^-^ In 



1 + 



2(ci +C2)fc 1 

'ci + k){c2 + k) e2*^"-l 



(34) 



As a check of this result, we note that the integral vanishes if ci = C2 — > oo or C\ — C2 — 0, 
thus reproducing the correct result for Dirichlet-Dirichlet and Neumann- Neumann boundary 
conditions. Dirichlet-Neumann boundary conditions (ci — > oo, C2 = 0) can also be treated 
exactly: in this case, the integral in Eq. (34) becomes 



I{a) := / 
Jo 



dkk'^ ^ Incoth/ca. 



(35) 



Integrating by parts yields 
2a /-"o „ k'^ 



dk 



sinh 2ka 



= i^E/;''"^»™- (2-^)^C(<i+l). (36) 



Inserting this result into Eq. (34) and using the identity V{2z) = (47r)~^/^ 2^^ r(^) r(^-M/2) 
[18] we finally obtain 



-d-l 



Eo(oo,0;a) 



1 

2^ 



(47r) 



-(d+l)/2 



C(rf+i), 



(37) 



which agrees with the correct result [17], thus giving us another check of Eq. (34). 

Next in simplicity are the following three cases: (i) C\ — C2 — c, (ii) ci = oo, C2 = c, and 
(iii) ci = 0, C2 = c. We shall denote them RR, DR and NR, respectively (R = Robin, D 
= Dirichlet, and N = Neumann). In all these cases, changing the variable of integration in 
Eq. (34) to q = k/c allows us to rewrite it as 



E^{c, a) = L^-^eSo.ica) {a = RR, DR, NR), 



(38) 
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where 



r(^)C(rf+i) 1 foo 



\^ fa{q) 



f,2qx _ I 



(39) 



with 



= ' -^Di^^^) = ' ■^NR(g) = ■ (40) 

The graphs of £a{ca) in three spatial dimensions are depicted in Fig. 1. We can conclude 
from it that the Casimir force between the plates: (i) is purely attractive in the RR case 
(i.e., Ci = C2 = c); (ii) is repulsive at short distances and attractive at long distances in the 
DR case, and (iii) is attractive at short distances and repulsive at long distances in the NR 

case. 

To understand those behaviors, let us consider a free field subject to Robin BC at 2; = 0, 
i.e., 920(x, 0) = c0(x, 0). If we write the 2;-dependent part of as ip{z) — sva.{kz + 5), the 
previous equation becomes tan 5 = k/c. It follows that 5 — > as /c — > 0, and 5 — > 7r/2 as 
k 00. In terms of ^p{z), this is equivalent to say that Robin BC tends to Dirichlet BC at 
low momentum, and to Neumann BC at high momentum. In the jargon of rcnormalization 
group theory, c = 00 is an infrared and c = is an ultraviolet attractive fixed point. 

Let us now return to Fig. 1. According to the analysis above, the RR curve should 
behave as the DD curve in the infrared (i.e., a — > 00) and as the NN curve in the ultraviolet 
(i.e., a — s> 0). In both cases, the Casimir force is purely attractive, and so it is in the RR 
case. The DR curve should behave as the DD curve as a ^ 00, and as the DN curve as 
a — > 0; indeed, this is what we observe: attraction at long distances and repulsion at short 
distances. The analysis of the NR curve is similar. 

Such considerations suggest an interesting possibility. Let us suppose that < Ci <S C2 < 
00. Then Eq{ci, C2; a) ~ Eq{oo, 00; a) as a ^ 00, and -£"0(^1, C2; a) ~ Eq{0, 0; a) as a — 0. In 
both these limits, therefore, the Casimir force is expected to be attractive. However, since 
the crossover from a Dirichlet-like to a Neumann-like BC takes place at different scales for 
each plate, there could be a range of distances for which one has a Dirichlet-like BC at one 
plate and a Neumann-like BC at the other, thus leading to a repulsive force between them. 
That such a possibility can indeed occur is shown in Fig. 2. 



IV. INTERACTING FIELD 

There is still another reason for the introduction of boundary terms in the action: the 
perturbative treatment of a renormalizable interacting theory requires boundary countert- 
erms in addition to the usual ones [8,12]. Thus, from a conceptual point of view, it is more 
natural to treat boundary conditions as resulting from the interaction of the field with a 
background than to impose them a priori. We shall illustrate this point with the calculation 
of the first order correction to the two-point Green function of the A0"^ theory in the presence 
of a flat boundary, where the field is submitted to Robin BC. 
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The renormalized Euclidean Lagrangian density of the theory is given by^ 

(^-\ (d.cPf + \ m'<j>' + cS{z) 0' + ^ 0' + >Cet, (41) 

with the field living in the half-space 2; > 0. For simplicity, we shall assume m — and 
c > 0. Cct contains the renormalization counterterms — the usual ones, 

SZ ,,9 Sm^ ,9 SX ,A 
>Cct,buik = -Y (9^^? + ^ + l! ' (42) 

plus boundary conterterms, which we shall exhibit later. 

As we have seen in Sec. 11, the unperturbed two-point Green function (Feynman propa- 
gator) is given by [see Eq. (18)] 

Go(x, x') := {ct>{x) 0(x'))o - / 7^ e^'^-^'^"^'^ ^o(^; ^, z% (43) 



where, in the massless theory. 



Golk; z, z') = ; h - smh( kz<) exp(-kz>). (44) 

c + k k 

If we neglect boundary counterterms, the first order correction to Qq is given by 

roo 

Qi{k]z,z') ^ - I dwgo{k;z,w)J2i{w)go{k;w,z'), (45) 
Jo 

where the one-loop self-energy, including the mass counterterm, is given by 

^^^^^ " ^ / (1^ ^ 

= / TTT^ + 7. + 5ml (46) 

(An ultraviolet cutoff A is imphcit in the integral above and wherever necessary below.) 
To fix Sm^ we impose that Tii{w) be finite for ■{/; > and vanishes for w — > 00 (the latter 
condition is equivalent to requiring that the physical mass is zero when one is infinitely far 
from the plate). These conditions imply 

^4/^1 (47) 

hence 



^The factor of \ in the boundary term in the action (4) is due to the fact that /o°° b{z) f{z) dz = 



9 



A f d'^q 



-2qw 



2 J {27rY \c + q 2q^ 



(48) 



Although is now a well defined function of w ior w > 0, this is not enough to 

ensure the finiteness of Qi. Indeed, ~ w''^^~^^ for ty — 0, causing the divergence of 

the integral in Eq. (45) for d > 2. This problem can be (partially) solved by the boundary 
counterterm 



4l!boundary = <^C5(Z) (f\ 



(49) 



It adds a new term to the self-energy Si(w), turning it into Si(w) = Si(w) + 26c d{w). 
Equation (45) then gets replaced by 

roo 

gi{k; z, z') ^ - dw go{k; z, w) Ei(u') Qo{k; w, z') - 5cQo{k] z, 0) ^o(^; 0, z'). (50) 
Jo 

Let us assume that Q < z < z'] then Eq. (50) yields 

Sdt, = -Sc _ _ / ^ (_ _ _j [/(O. .) + /) + /(.'. oo)]. (51) 

where 

I{a,b):^ dwe-''i'^go{k;z,w)go{k;w,z'). (52) 

Ja 

Since we are interested in the large-g behaviour of the integrand, we may restrict our atten- 
tion to /(O, z) — I{z, z') and I{z', oo) are exponentially supressed in that limit. Performing 
the integration, we obtain 



l{0,z) 



-k{z+z' 



^0 



dw e 



-2qw 



q— >oo 



(c+fc)2 



1 c 
2^ + 2? + ^ 



+ — smh(kw) 

c+k k ^ ' 

' 1 ^ 



(53) 



If we further multiply this result by the term in parentheses in Eq. (51), we conclude that 
the integrand in that equation behaves for large q as 



-fe(z + 2') 



{c^kf 
Therefore, if we choose 



q— >oo 



{c + ky 



4g(g + c) 



or 



(54) 



d'^q 



8 J {27iY q{q + c) 



(55) 



with 5c finite, Qi{k\ z^ z') becomes finite for < 4. 

To see why the boundary conterterm (49) only solves the problem partially, let us take 
2; = 2;' = in Eq. (50). A straightforward calculation then shows that 
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1 



Sc 



Sc + 



A f d'^q 



1 



1 



1 



2 J {2t:Y \c + q 2qJ 2{q + k) 

\ r d<^q 1 



c + A;)2 8{c+k) J {2nY q{q + c){q + k)' 



(56) 



The result is finite for d = 2, but diverges for c? = 3. In this case, another boundary 
counterterm is needed. As we show below, it is given by 



4t?boundary = ^bS{z) (t){dn - c)4), 



(57) 



with dn the interior normal derivative (in our case, dn = d^)- Indeed, such a counterterm 
gives rise to an extra contribution to Qi{k; 0, 0), given by 

Agi{k; 0, 0) = -5b / dw 5{w) {d^ - 2c) g^{k; 0, w) 
Jo 



= -Sb 
5b 



roo 

/ dw 5(w) (dw — 2c) 
Jo 



-2kw 



[c + kf 



c + k 



If we choose 



5b 



8 J 



A f d'^q 



+ 5b, 



{2n)'^q^{q + c) 

with 5b finite, we cancel the divergence in Qi{k] 0, 0) for c? = 3. Indeed, 



(58) 



(59) 



^;i(A;;0,0) + A6;i(A;;0,0) 



c + A; 
1 

c + k 



5b- 
5b- 



5c 
c + k 
5c 



Xk r d'^q 
~Y J (2^ 



{2n)d q2(q + c){q + k) 



+ 



A 



C + k 167r^ c — k \c 



{d = 3). (60) 



The boundary counterterm (57) is ineffective if z,z' > 0; in this case, A^i(/c; z') is 
identically zero. (As a consequence, Qi{k\z,z') is insentive to the choice of 5b if z^z' > 0.) 
On the other hand, that counterterm is necessary if z = and z' > (or vice- versa); in 
this case, it is possible to show that the same choice (59) for 5b also ensures the finiteness 
of gi{k]z,z'). 

In order to fix the finite part of the boundary counterterms, 5b and 5c, one has to impose 
a pair of renormalization conditions. A natural choice, for it is satisfied at tree level, is given 

by 



g{K;0,0) = {c + k)-\ 



(61) 



^^(.;0,0) 



Ic + k) 



(62) 



where k is some arbitrary, but nonzero, mass scale. 
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We end this section with some remarks: 

(A) Because of the boundary counterterms, the renormahzed two-point Green function 
does not satisfy Robin BC at ^ = — except for c — > oo (Dirichlet BC); in this case, the 
boundary condition is preserved at each order in perturbation theory [12]. 

(B) Inclusion of another plate at 2; = a does not affect significantly the overall picture. 
New ultraviolet divergences arise, but the theory is made finite with the same type of 
counterterms used in the case of a single plate. 

V. CONCLUSIONS 

In this work we have computed the Casimir energy of a free massless scalar field subject 
to independent Robin boundary conditions on two parallel plates in d spatial dimensions. 
It was shown that for mixed Dirichlet-Robin BC the Casimir energy as a function of the 
distance a between the plates displays a minimum. We managed to understand the behavior 
of the Casimir energy as a function of c relying on an analogy with the renormalization 
group flows expected from infrared/ultraviolet fixed points. It was found that the Dirichlet 
BC (c — > 00) is analogous to an attractive infrared fixed point whereas the Neumann BC 
(c = 0) resembles an attractive ultraviolet fixed point. This interpretation is consistent with 
the numerical results shown in Fig. 1 for the Casimir energy as a function of ca, and suggests 
that a crossover from Dirichlet-like to Neumann-like behaviors at the plates may lead to a 
repulsive force between them. 

We also provided a detailed analysis of renormalization for the two-point Green function 
x') at first order in A in the theory. For simplicity we worked out the case of 
Robin BC at a single flat boundary. We have shown that, in addition to the usual "bulk" 
counterterms, one boundary counterterm is necessary to render G finite in the bulk, and a 
second one is necessary if at least one of its arguments lies on the boundary. This analysis 
is a necessary step in the computation of radiative corrections to the Casimir energy, which 
we intend to present elsewhere. 

In view of the widespread interest on the Casimir effect and its possible technological 
applications it is worthwhile to seek alternative computational tools which may go beyond 
or complement the existing ones. For instance, in [19] a resummation scheme to compute 
the leading radiative corrections to the Casimir energy was suggested. On the other hand, 
the method outlined in the present work seems well suited to the computation of radiative 
corrections either in a perturbative setting or eventually via application of semi-classical 
methods. Its starting point is the indirect implementation of the BC by means of appropriate 
terms in the action functional, which is then reexpressed in terms of two kinds of fields: a 
field ??(x, satisfying Dirichlet BC at the surfaces, and two surface fields 0j(x) localized 
on the planes, depending only on the remaining d transverse coordinates. The advantage 
of this procedure is that the functional integration over the surface fields is unconstrained, 
i.e., one does not have to enforce explicitly the Robin BC on the fields. 

The implementation of BC via local terms in the action is usually employed in studies 
of boundary critical phenomena [8]. In that context, it can be shown that Dirichlet and 
Neumann BC correspond to the so-called ordinary (c — > 00) and special transitions (c = 
0), respectively. The Robin BC is relevant in the study of the crossover between those 
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universality classes, for which, however, the computations become much more involved. It is 
relevant also for the analysis of the ordinary transition; in this case, however, one may resort 
to an expansion in powers of [8]. We expect that the methods proposed here can be 
useful to the study of the crossover for the relevant case of two fiat planes. This is presently 
under investigation. 
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FIGURES 
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FIG. 1. £a{ca) [see Eq. (38)] in three spatial dimensions for a = RR {ci = C2 = c), a = DR 
(ci = oo, C2 = c), and a = NR (ci = 0, C2 = c). 
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FIG. 2. Casimir energy per unit area vs. distance between plates (ci = 1/2, C2 = 2, d = 3). 
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